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THE  MIXINa  OP  JETS  IN  A  CHANNEL  WITH 


VARIABLE  CROSS-SECTION 
0.  V.  Yakovlevs kiy 


In  gas  machines  we  encoxmter  currents  with  non- 
uniform  profiles  of  geut-dynamlc  parameters,  e.g. , 
velocity,  this  non-unlformlty  gradually  smoothing 
out  downstream.  Similar  flows  are  observed  In 
pa3?tlovilar  In  the  mixing  chambers  of  ejectors,  two- 
clrcult  Jet  engines  [1],  etc.  In  some  practical 
cases  these  flows  occ\ir  In  channels,  the  form  of 
which  differs  from  the  cylindrical,  l.e..  In  chan¬ 
nels  with  variable  cross-section. 

In  planning  such  gas  machines  and  apparatus  It 
Is  necessary  to  know  the  laws  of  change  In  the  non- 
unlformlty  of  the  stream  along  the  flow-through 
section.  The  methods  of  calculating  the  parameters 
of  a  non-uniform  stream  In  a  channel  of  arbitrary 
shape  have  not,  however,  yet  been  established. 

There  exists  only  the  monograph  [2]  of  A.  Ya.  Cherkez 
on  a  theory  devoted  to  flow  In  the  cylindrical  blend¬ 
ing  chamber  of  an  ejector. 

In  the  present  work  an  attenqpt  Is  made  to  devise  a 
method  of  calculating  the  flow  of  a  fluid  with  non- 
unlfom  velocity  distribution  of  the  Jet  type  In  a 
channel  with  variable  cross-section. 


1.  Initial  premises  and  equatlcms.  When  mixing  coaxial  streams 
In  a  channel  (Pig.  1)  It  Is,  as  Is  well  known  [2],  possible  to  sepaiv 
ate  at  least  two  essentially  different  regions  of  flow.  In  the  first 
of  these,  called  the  primary  zone,  the  dimensions  of  the  region  of 
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flow  mixture  are  less  than  the  transverse  dimensions  of  the  channel; 
whereas  in  the  following  secondary  zone  the  mixing  process  already 
embraces  the  whole  cross-section  of  the  stream.  The  flow  in  pre¬ 
cisely  this  zone  will  be  examined  below. 

From  analysis  of  the  velocity  profiles  in  the  second  zone  of 
the  stream  in  cylindrical  and  conical  (ef fuser)  channels  the  con¬ 
clusion  has  been  successfully  drawn  that  the  the  velocity  distribu¬ 
tions  in  different  cross-sections  reduced  to  dimensionless  form  by 
two  parameters  maintain  their  similarity  along  the  channel.  Con¬ 
cretely,  it  is  question  of  autos Imulat ion  of  the  profile  of  dimen¬ 
sionless  velocity  of  the  form 

..•=^,.^(1)  (1.1) 

Here  u  is  the  longltudliwil  component  of  flow  velocity;  veloc¬ 
ity  on  the  X  axis;  £,  average  (areal)  flow  velocity  in  the  given 
section;  Rj  present  radius  of  the  channel;  and  2,  distance  from  the 
axis. 

The  e3q>erlmental  data  of  A.  Ya.  Cherkez  [2]  and  K.  Vllctorin  ['^] 
for  a  cylindrical  channel  are  shown  in  Fig.  2  treated  as  the  function 
where  |  ■  j/R}  and  the  results  of  the  author’s  experiments  for 
a  mixing  chamber  of  the  ef fuser  type  (truncated  cone  with  aperture 
angle  of  8^20*),  in  Fig.  3.  As  we  see,  the  profile  of  dimensionless 
velocity  du^(0  is  similar  in  both  cases  in  the  different  cross-sec¬ 
tions  of  the  stream. 

Fig.  i.  Piagram  of 
Jet  propagation  in  a 
channel  with  variable 
cross-section. 
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Pig.  2.  Profile  of  relative 
excess  velocity  6u®(0  In  a 
cylindrical  channel.  Ctirve 
from  Ponmila  (l.2]j  points 
from  ej^erlments  (dark^lr 
ejector  [2],  llght-icater 
mixer  [3])  for  following  dis¬ 
tances  from  entry  cross-sec¬ 
tion  relative  to  channel 
radius:  l)  7.5*  2) 

3)  13.5,  4)  7.3,  5) 

6)  9.4,  7)  10.4. 


1*5, 

3, 


Pig.  3.  Profile  of  rela¬ 
tive  excess  velocity  6u®(4) 
In  effuser  channel  (a  - 
=  -0.0721).  Cxirve  from 
Pormula  (1.2);  points  from 
author's  experiments  at  fol¬ 
lowing  distances  from  entry 
cross-section:  1)  462  mm, 

2)  662  mm,  3)  746  mm. 


By  analogy  with  the  free  turbulent  stream  let  us  assiune  that 
autoslmulatlng  distribution  of  dimensionless  velocity  In  the  cross- 
section  of  the  Jet  stream  In  the  channel  Is  described  by  the  so-called 
Schllchtlng  law  which  Is  written 


(1.2) 


when  applied  to  the  case  \znder  examination. 

Coefficient  A(l)  Is  a  numerical  constant  which  Is  selected  so 
that  the  equation  of  discontinuity  will  be  Identically  satisfied  when 
using  profile  (1.2). 

The  ctirves  corresponding  to  relationship  (1*2)  are  also  shown  In 
Pigs.  2  and  3;  It  Is  evident  that  In  their  first  approximation  they 
describe  the  autoslmulatlng  profile  of  velocity  du^(0  rather  well. 

Thus,  one  of  the  Initial  prerequisites  of  the  calculation  method 
Is  the  similarity  In  the  profiles  of  dimensionless  velocity  which  was 
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discovered  on  analyzing  the  experimental  data;  In  this,  the  autoslmu- 
latlng  distribution  of  velocity  6u®(0  may  be  described  by  relation¬ 
ship  (1.2). 

Moreover,  the  folloirlng  assumptions  were  made  In  the  calcula¬ 
tions:  the  fluid  Is  considered  Inconpresslble  (p  »  const),  the  static 
pi?essui>e  £  assumed  to  be  constant  at  every  cross-section,  and  the 
force  of  friction  between  the  fluid  and  the  wall  of  the  channel  Is 
neglected.  The  last  a8suDQ>tlon  Is  Justified  In  that  the  extent  of 
the  channel  In  which  mixture  takes  place  Is  usually  small  and  amounts 
to  only  a  few  diameters  of  the  channel. 

Consequently  the  problem  Is  reduced  to  determining  the  dependence 
of  three  tuiknown  values,  u^,^,  g,  and  p,  on  longitudinal  coordinate  x. 
Axlsymmetrlcal  fluid  flow  In  the  examined  case  Is,  as  Is  well  known 
[4],  determi^d  by  the  equations 


(1.3) 


Here  u  and  v  are  the  o<»q}onents  of  velocity  of  average  motion; 
and  the  stress  of  turbulent  friction  according  to  L.  Prandtl  Is 
described  by  the  fozmula 


Tt 


System  (1.?)  may  be  easily  omiverted  to  the  form 

If  both  halves  of  the  first  of  the  aquations  (1.4)  are  multiplied 
by  dy  axid  Integrated  from  0  to  an  arbitrary  and  then  the  valiae  pvy 
from  the  second  of  these  aquations  (1.4)  is  substituted,  we  get 

•  t 

This  is  the  first  of  the  correlations  of  the  system  of  equations 


for  detezmlnlng  the  above-listed  three  variables.  The  second 


coxrelatlon  Is  derived  from  (1.5)  If  the  upper  limit  In  the  Integrals 
entering  here  Is  set  equal  to  channel  radius  R  and  the  boundary  con¬ 
dition  at  the  wall,  y  »  R,  ■*  0,  Is  satisfied: 


(1.6) 


The  third  correlation  In  our  system  Is  the  condition  for  main¬ 


taining  mass  In  the  channel 


A 

gj^pHydy  =  0,  Or  C=r  const  (1.7) 

0 

and  serves  to  determine  tte  average  (areal)  stream  velocity  In  each 
cross-section. 

Let  vus  Introduce  the  dimensionless  coordinates  and  dimensionless 
velocities 


_  *  e  _  V 
R  '  »  — 


Jf": 


A  u  a 


=  (1-8) 


where  R  Is  the  radlvus  of  the  channel  at  some  fljced  cross-section. 

Let  us  replace  the  gradient  of  pressure  In  Eq.  (1.5)  by  Its 
expression  In  (1.6);  In  doing  so  we  will  assume  that  the  relationship 
of  bilxture  path  1  to  the  channel  raullus  In  the  same  section  remains 
constant  along  the  channel  and  can  depend  only  on  l.e.,  1/R  >  x(Ol 
In  addition,  we  will  use  expression  (1.2)  to  distribute  velocity  u 
In  the  cross-section  of  the  stream. 

As  a  result  we  obtain  the  following  differential  eqtiatlon  for 


determining  the  dimensionless  velocity  on  ohanxiel  axis  U 


(1.9) 


Here  the  top  sign  before  the  first  member  on  the  right  side  Is 
taken  If  the  velocity  on  the  flow  axis  Is  greater  than  at  the  wall; 
and  the  bottcai  sign.  In  the  opposite  case.  Moreover,  In  Eq.  (1.9) 
the  following  designations  are  Introduced: 
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^(6).-2[fia)-6*i?(i))-ji/a)i3.i(i)+9(6)i  (1.10) 

pa)^Ba)-VB(i)-2A(i)M(i) 

A{l)~2\if(l)ldl,  il(l) -0.2571,.  .B(l)-0.1335 

0  .  • 

If  function  q>(0  were  an  autoslnulatlng  solution  of  the  motion 

equations  Eq.  (1.9)  would  not  contain  that  Is,  the  coefficients 

of  M(C),  N(0«  P(^)  would  be  proportional  to  each  other. 

But  In  our  case  these  coefficients  possess  no  exact  proportionality; 

one  can  convince  himself  of  this  from  an  examination  of  Fig.  4  which 

lists  the  relatlcmshlps  M((),  N(0«  and  P(0*  ^  this  we  do  not 

examine  the  change  In  valiie  of  T(0  since  it  has  in  principle  no 

2 

significance  as  it  is  the  co-factor  of  the  empirical  constant  h  .  In 
order  to  simplify  the  subsequent  analysis  we  will  relate  the  coeffl«» 
dents  of  Eq.  (1.9)  to  the  value  of  M(0  and  Introduce  the  following 
designations  for  these  relative  coefficients; 


Substituting  these  coefficients  In  Eq.  (1.9)  and  also  introduc¬ 
ing  the  val\ie  k  ■  hS)(0  we  finally  get 

(1.11) 


IP(l+aU)^^ZfkU*  +  2^U(2  +  eU) 

Since  there  are  no  objective  criteria  for  the  selection  of  the 

value  of  $  at  which  the  coefficients  of  a(()«  I’CO  and  c(0  should 

be  figured,  we  will,  in  our  desire  to  avoid  arbitrariness,  find  their 

average  integral  values  In  accordance  with  the  formulas 

1 


rf/J* 


Cosipleting  the  integration  we  finally  get 

0.34W,  15.0, 

c^-  0.1664 


(1.12) 
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We  will  use  these  values  of  the  coefficients  of  Eq.  (l«9)  In 
what  follows. 

2.  Analysis  of  the  Solution.  In  the  general  case  when  the 
relationship  R®(x°)  has  an  arbitrary  form,  Eq.  (l.ll)  Is  an  Abelian 
equation  of  the  second  sort  [3]  and  Its  solution  Is  not  expressed  In 
quadratures.  In  two  frequent  cases  which  are  examined  below^  the 
solution  of  differential  Eq.  (l.ll)  may  be  successfully  expressed 
through  elementary  functions. 

Case  1.  A  cylindrical  channel  (R^  »  const  =  l).  In  the  case 
of  a  cylindrical  channel  Eq.  (l.ll)  Is  greatly  simplified  and  Its 
solution  In  the  boundary  condition  where  x°  =  0,  U  =  U*  has  the  form 

If  we  substitute  here  the  numerical  value  of  coefficient  a 
according  to  Fozmila  (1.12)  and  relate  the  excess  velocity  on  the 
stream  axis  U  to  Its  value  In  the  first  section  l.e..  Introduce 
U®  =  V/U^,  we  finally  get 

^-l-0.80U,lgU^^±kUX  (2.1) 

The  eiiq;)lrleal  constant  k  In  this  equation  Is  approximately  0.26 
according  to  estimates  based  on  the  experimental  data  of  different 
authors  [2,  3,  6].  The  change  In  the  static  pressure  along  the 
channel  Is  found  from  Eq.  (1.6).  The  change  In  static  pressure  Is 
related  to  velocity  pressure  corresponding  to  the  average  (areal) 
velocity  £,  which  In  the  case  examined  remains  constant  along  the 
channel  (g  -  const  •  and  Is  determined  from  the  formula 

A/,*  =  -0.0674  (2.2) 

This  fozmila  enables  us  to  find  the  distribution  Ap®(x®)  since 
the  relationship  n®(x®)  Is  known  and  Is  detezmlned  by  Eq.  (2.1). 
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The  change  In  the  value  of  In  a  cylindrical  chanxiel  with  a 
different  pattern*  of  velocity  profile  In  the  initial  section  charac¬ 
teristic  of  excess  velocity  Is  Illustrated  In  Fig.  3-  It  Is 
evident  that  the  greater  the  Initial  nonunlfonolty,  the  more  Intense 
Is  Its  equalization  also,  dtirlng  which  In  accordance  with  Fonoula  (2.2) 
a  greater  Increase  In  presstire  also  occurs. 

Case  2.  A  channel  with  vari¬ 
able  cross-section  and  a  straight- 
line  generatrix  (R®  *  1  +  ox®). 

The  differential  Eq.  (i.ll)  In  the 
case  In  question  takes  the  form 

(1  +  ax^)  (1  +  aU)  =  AaU  +  (2M^k)U* 


/ 

\ 

— 

‘''^4 

\ 

Fig.  4.  Change  In  the  coef¬ 
ficients  M(0#  and  P(0 
In  cross-section  of  stream. 


o 
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The  solution  of  this  equation  with  the  boundary  condition  of 
0  and  U  -  Is  expressed  by  the  following  single  formula: 


(2.5) 


Ifere  ajra  Introduced  the  following  designations 


m 


If,  as  In  the  preceding  case,  we  Introduce  the  relative  axial 
velocity  U®  relationship  (2.5)  will  be  rewritten  thus: 


*  The  nonunlfoxTolty  of  the  velocity  profile  may  be  characterized 
by  the  coefficient  v  -  1  -t,  where  the  value  of  t  Is  the  ratio  of  the 
velocity  averaged  for  flow  rate  to  the  velocity  averaged  for  area 

R  R 

T-2a»Jii»|frfr/(2ji»wr)* 

Ad(^ting  Foxnula  (1.2)  for  the  distribution  of  the  velocity  we 

can  show  that  t  ■>  1  4-  O.O674  U^.  Thus  the  coefficient  of  nonunlformlty 

V  -  0.0674  proves  to  be  Immediately  connected  with  the  excess  axial 
velocity  TJ. 
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(2.4) 


i 


As  ire  see,  the  change  In  velocity  al(»ig  the  axis  of  the  channel 
Is  determined  hy  two  parameters,  and  a,  characterizing  the  Initial 
pattern  of  the  velocity  profile  and  the  gemnetrlc  properties  of  the 
channel,  respectively. 

Knowing  the  relationship  IJ^(x^)  we  may  find  the  distribution  of 
the  static  pressure  along  the  channel  from  £q.  (1.6).  Omitting  the 
Intermediate  calculations  we  will  write  the  final  es^resslon  relating 
the  change  In  pressure,  as  In  the  preceding  case,  to  velocity  pressure 

O 

pg  /2  corresponding  to  the  average  (areal)  stream  velocity  In  cross- 
section  x^  e  0. 

Here 

0.1348  (/,  ^  0.0674  (3  —  4«)  (f,* 

*-m{l-4n)*  ‘ 

The  evolutlm  of  the  nonuniformity  of  the  velocity  profile  In 
channels  of  the  effuser  and  diffuser  type  Is  shown  In  Fig.  3  In  the 
form  of  the  relationship  U®(x®)  for  three  values  (40.05,  -0.05,  -0.1) 
for  parameter  a  corresponding  to  central  angles  of  cones  of  5^40* 
and  11^50*.  As  should  have  been  e3q>ected  the  most  rapid  damping  of 
the  nonuniformity  occurs  In  the  effuser.  In  which  process  the  effect 
of  parameter  becomes  weaker  with  the  Increase  In  the  degree  of 
constriction  of  the  channel,  and,  for  exaaple  when  a  -  -0.10  the 
curves  of  II^(x^)  for  «  0.5  and  1.0  practically  coalesce  Into  one. 

The  change  In  the  nonunlfozmlty  of  the  velocity  profile  along 
the  channel  occurs  under  the  Influence  of  two  factors:  turbulent 
mixing  (proportional  to  the  cross-sectional  gradient  of  the  velocity) 
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always  leading  to  the  leveling  of  the  nonuniformity,  and  unequal 
acceleration  of  the  fluid  In  different  flow  tubes  on  changing  the 
area  of  the  cross-section  (with  the  same  pressure  change  the  slower 
Jets  accelerate  In  the  effuser  or  slow  down  In  the  diffuser  more  than 
do  the  faster  Jets).  The  first  of  these  factors  In  Its  p\irest  form 
appears  In  the  exanqiles  examined  In  Fig.  3  In  the  case  of  a  cylindri¬ 
cal  channel  (a  -  0),  and  the  secwid  factor  begins  to  prevail  most 
powerfully  over  the  first  In  the  rapid  area  change  of  the  cross- 
section,  that  Is,  In  a  large  gradient  of  pressure  (this  Is  observed 
In  the  exaiiq>les  In  Fig.  5»  e.g. ,  for  the  effuser  even  when  a  >  -O.iO). 
In  the  flow  In  the  effuser  both  factors  act  In  the  same  direction  and 
In  the  flow  In  the  diffuser  the  second  factor.  In  distinction  from 
the  first,  favors  the  Increase  of  velocity-profile  nonuniformity.  The 
common  Inflxienoe  of  the  noted  factors  Is  very  clear  In  the  cases  of 
flow  shown  In  Fig.  3  In  a  weak  diffuser  (central  angle  3^40*)*  At  a 
greater  Initial  nonuniformity  (U^  -  l.O)  both  factors  are  approximately 
the  same  In  stwngth:  the  ciirve  of  U®(x®)  diverges  weakly  from  the 
line  of  »  1  (In  the  direction  of  U  <  l);  but  at  a  smaller  non- 
unlformlty  (U^  *0.3)  the  first  of  the  mentioned  factors  becomes 
considerably  weaker  than  the  second  and  as  a  resiilt  of  the  unfavorable 
effect  of  the  second  factor  the  nonMnlfozmlty  of  the  velocity  profile 
proves  to  be  pn^p^esslve  downstream  (U®  >1). 

In  flow  calculations  In  conical  diffuser  channels  we  should 
keep  In  mind  the  possibility  of  the  stream  breaking  away  from  the 
wall  of  the  channel  when  the  size  of  the  central  angle  exceeds  8  to 
10®.  The  ealoulatl<m  method  proposed  here  In  the  case  of  break-away 
diffusers  may  be  applied  only  In  the  flow  region  where  the  zone  of 
reverse  oiunents  has  already  disappeared  and  the  velocity  profile  of 
fimn  (l.l)  has  beoosm  autMlmulatlng. 
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3«  Comparison,  with  the  experimental  data.  Research  on  flow 
with  a  Jet  profile  of  velocity  In  a  channel  of  cylindrical  shape  has 
been  conducted  hy  several  authors.  In  particular,  data  of  Interest 
to  us  are  contained  In  works  hy  A.  Ya.  Cherkez  [2]  and  also  hy 
W.  Richards  and  V.  Oshome  [6]  on  low-pressure  air  ejectors  and  In 
reseatpch  hy  V.  M.  Papin  on  a  water  ejector  [7].  Figure  6  shows  the 
experimental  res\xlts  of  these  authors  together  with  the  TJ®(x®)  curve 
cooputed  of  the  theoretical  relationship  (2.1).  In  doing  this.  In 
each  case  we  selected  as  a  heglimlng  for  the  calculation  on  the  x-axls 
the  cross-section  In  which  =  1;  we  took  the  value  of  the  esplrloal 
coefficient  k  to  he  0.26.  As  we  see,  the  calculation  method  proposed 
here  leads  to  conpletely  satisfactory  agreement  with  the  experimental 
data  on  hand.  The  constancy  of  the  enplrlcal  coefficient  In  the 
calculation  formula  under  different  experimental  conditions  which  was 


discovered  from  a  conpauclson  of  the  theory  with  e3q)erlment  should  he 


especially  enphaslzed. 


Fig.  3*  Calculated  rela- 
tlmuihlp  iD®(x®)  for  dif¬ 
fuser  (a  ■  0.05),  cyllj>> 
drloal  tube  (a  >*  O),  and 
effuser  channel  (a  - 
-  -0.05,  -0.10).  Solid 
curves  for  U*  -  0.5; 

dashed,  for  -  i.O. 
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6. 

Conparlson  of  ci 

otilated  curve  of  dancing  of 
relative  axial  velocity 
with  e:q)erlmental  data  of 
various  authors.  Curve  from 
Foxvula  (2.1);  points  from 
experiments  of  1)  Cherkez 
[2J,  2  and  3)  Richards  and 
Oshome  [61,  <;)  Papin  [7]. 
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Fig.  7.  Change  In  axial 
velocity  of  non- 

uniform  stream  In  effuser 
with  Btralght-llne  gener¬ 
atrix  (a  -  -0.0721). 

Curve  from  calculation, 
points  from  author’s 
esiperlments. 


Fig.  8.  Change  In  static 
press\ire  Ap®(x®)  In  non- 
uniform  stream  In  effuser 
with  Btralght-llne  gener¬ 
atrix  (a  -  -O.O72I). 

Ciurve  from  calcvilatlon, 
points  from  author's 
experiments. 


It  was  of  great  Interest  to  c^i^are  the  theoretical  and  esqperl- 
mental  results  In  the  case  of  a  channel  with  variable  cross-section. 
Such  experlmeiits  were  conducted  by  the  author  on  a  model  of  a  effuser 
mixing  chamber  idilch  represented  a  tnincated  cone  with  an  angle  of 
8^20'  (a  »  -0.0721)  1040  mm  long  and  with  a  cross-sectional  diameter 
of  100  mm.  Two  coaxial  streams  of  air,  the  flow  rates  per  second  of 
i^eh  could  be  Independently  regulated,  were  supplied  to  the  entx>y 
section  of  the  mixing  chamber;  this  allowed  us  to  create  different 
nonunlfonnltles  of  velocity  profile  In  the  entry  section  of  the 
channel.  Figure  7  shows  the  data  of  one  of  these  experiments  In  the 
foxm  of  relatlmishlp  U^(x^};  the  calculation  curve  from  Fonnula  (2.4) 
Is  also  shown  here.  The  noniinlformlty  of  the  stream  In  the  section 
taken  to  be  Initial  was  characterized  by  a  value  of  «  1.0.  The 
«qplrloal  coefficient  k  was  taken  In  the  case  examined  to  be  0.2$, 
l.e.,  equal  to  the  value  derived  for  the  flow  In  a  cylindrical  tube. 

Figure  8  shows  the  e:q)erlmental  data  obtained  In  the  author's 
experiments  on  the  distribution  of  static  pressure  along  the  wall  of 
the  channel;  here  Is  also  entered  the  curve  defined  by  relatlonshlpB 
(2.4)  and  (2.3)*  From  an  analysis  of  Figs.  7  and  8  It  follows  that 
the  caloulatlm  method  pressed  here  Is  completely  suitable  for 
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practical  uae  In  studying  the  motion  of  a  fluid  with  a  nonunlform 
velocity  pi?oflle  of  the  Jet  type  In  a  effuser  channel  with  a  straight- 
line  generatrix. 


KKJi'JbihENCT^-*’! 


g«  N.  Abramovich.  Applied  Oas  Dynamics.  Oostekhlzdat,  1955. 
2.  0.  N.  Alaramovlch.  Theory  of  Turbulent  Streams.  Plzmatglz, 


5*  yhtersuchung  turbulenter  Mlschvorg&nge.  Forsch. 

Ing.-wes.,  Hd.  12.  s.  16-30,  1941. 

^  Qol|<lflhteyn.  ed.  The  Present  State  of  the  Hydro-Aero¬ 

dynamics  or  a  VISCOUS  Fluid,  Vol.  1.  IL,  1948.  [Translation  of 
Modem  Developments  In  Fluid  Dynamics,  Vols.  1  and  2.  Oxford  Uni¬ 
versity  Press,  1938]. 


5*  E.  Knmirw.  Handbook  of  Ordinary  Differential  Equations. 
Flzmatglz,  1901. 


tor. 


W.  Q.  Richards,  W.  C.  Osborne, 
mst.  Heat  and  vent liat.  Engrs 


The  very  low  pressure  ejec- 
.,  Vol.  27,  pp.  172-179,  1959. 


7»  V.  M.  Panin. 
Occurring  In  the  Mixing 
84,  No.  5,  1952. 


Experimental  Investigation  of  Phenomena 
Chamber  of  a  Water-Jet  Piaap.  DAN  SSSR,  Vol 


Received  April  I3,  196I. 


FTD-Tr-62-15ri-l+2+4 


-13- 


DISTRIBUTION  LIST 


) 


DBPAR'B'iENT  OP  LSFENSE  Ifr.  Copies  I£AJ03  AIR  CCBilklANDS  Nr.  Copies 


APSC 

SCPTR 

ASTIA 

HBADQU^IRTERS  USAP 

TD-Bla 

APCi;f-3D2 

1 

TD-Blb 

AEDC  (AEY) 

AilL  (ARB) 

1 

SSD  (SSP) 
BSD  (BSP) 
APPTC  (PTY) 
APSWC  (SWP) 

OTHER  AGENCIES 

1 

25 

5 

3 

1 

2 

1 

1 

1 


CIA  1 

USA  6 

AID  .  2 

OTS  2 

ABC  2 

p;/s  1 

NASA  1 

RAND  1 


V 


p:d-tt-  62-1571/1+244 


14 


